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SUMMARY 


The problem of sidewall boundary-layer effects in airfoil testing is 
treated by considering the changes in the flow area due to boundary-layer 
thinning under the influence of the airfoil flowfield. Using von Karman's 
momentum integral equation, it is shown that the sidewall boundary-layer 
thickness in the region of the airfoil can reduce to about half the undis- 
turbed value under the conditions prevailing in testing of supercritical 
airfoils. A Mach number correction due to this increased width of the flow 
passage is proposed. Using the small disturbance approximation, the effect 
of the sidewall boundary-layers is shown to be equivalent to a change in the 
test Mach number and also in the airfoil thickness. Comparison of the re- 
sults of this approach with other similarity rules and correlation of the 
experimental data demonstrate the applicability of the analysis presented 
from low speeds to transonic speeds. 

INTRODUCTION 

Numerous studies have been reported in the literature to assess the 
wall interference effects in wind tunnel testing of models. These studies 
are generally based on inviscid calculations and the influence of wall 
boundary-layers on interference effects often ignored. However, in testing 
of airfoil sections in two-dimensional wind tunnels, recent studies at ONERA 
(ref. 1) have indicated that the presence of the boundary layers on the 
sidewalls of the tunnel can influence the measurements even on the airfoil 
midspan. To keep this influence within acceptable limits, it is desirable 
to use suitable sidewall boundary layer control. In the absence of such a 
control, appropriate corrections are to be made to the measured data to 
account for the effects of the sidewall boundary layer. Particularly, with 



the current development of flexible walls to nearly eliminate top and bottom 
interference in airfoil testing, the necessity to assess the extent of side- 
wall boundary-layer influence is a matter of prime concern. 

Methods to account for the sidewall boundary-1 ayer effects have been 
primarily based on two different considerations. Firstly, as proposed by 
Preston (ref. 2) in 1944, it is assumed that the trailing vorticity due to 
loss of lift in the boundary-1 ayer is shed and this in turn causes a change 
in the effective angle of incidence. In the second approach, proposed inde- 
pendently by Barnwell (References 3 and 4), and Winter and Smith (Reference 
5), the effects are attributed to changes in the sidewall boundary-layer 
thickness due to the model induced pressure field. In the Barnwell's 
method, based on the small disturbance approximation both the sidewall 
boundary-layer effects and the compressibility effects are combined into a 
single factor similar to Prandtl-Gl auret rule. This was later extended to 
transonic speeds by Sewall (Reference 6) by using the von Karman similarity 
rule. 

In the Barnwel 1 -Sewall method, it is assumed that the sidewall bound- 
ary-layers induce a spanwise velocity which varies linearly across the test 
section width. Recent investigations (Reference 7) on a super-critical 
airfoil in the NASA Langley 0.3-m Transonic Cryogenic Tunnel have indicated 
that correction to the test Mach number based on the Barnwel 1-Sewal 1 method 
accounted for most of the blockage effects and was sufficient to give a good 
match between the measured pressure distribution and the predictions of an 
advanced computer code. Particul arly, on supercritical airfoils, the pres- 
sure distribution is nearly constant over the forward portion and if the 
growth of the sidewall boundary-1 ayer due to viscous effects is small, it is 
reasonable to expect the sidewall boundary-layer thickness in this region to 
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be nearly constant at a value lower than the undisturbed thickness. Hence 
the primary effect of the sidewall boundary-layers appears to be due to 
change in the area of the flow passage rather than the presence of signifi- 
cant spanwise velocities across the width of the tunnel. In the present 
report, a new method to account for the change in the effective free-stream 
Mach number due to change in the sidewall boundary-layer thickness is pro- 
posed. Further, it is demonstrated using the small disturbance approxima- 
tions, that the effect of the sidewall boundary-layers can be considered as 
causing changes in both the free stream Mach number and the thickness of the 
airfoil . 
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NOMENCLATURE 

Area 

Speed of sound 

Width of the two-dimensional tunnel 
Airfoil chord 
Pressure coefficient 

Shape factor of the sidewall boundary-1 ayer 

* 

i 25 

Constant, = (2 + — - M 2 ) — — 

H °° b 

Local Mach number 
Corrected Mach number 

Effective Mach number in the equivalent 2-D flow 
Free-stream Mach nimber 

Correction for freestream Mach number (= M - M ) 

c 00 

Local velocity 

Effective free-stream velocity in the equivalent 2-D flow 
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U, Free-stream velocity 

u Perturbation velocity 

Perturbation velocity due to sidewall boundary-layers 

u Perturbation velocity in the wind tunnel 

w 

x, y Streamwise and normal coordinates 

Y Ratio of specific heats 

* 

6 Boundary-layer displacement thickness 

★ 

6 Undisturbed value of displacement thickness 

0 Boundary-1 ayer momentum thickness 

S , n Transformed coordinates 
p Density 

t Airfoil thickness ratio (geometrical) 

Effective thickness in the equivalent 2-D flow 

t Wall shear stress 

w 

<j> Velocity potential 

<t» w Velocity potential corresponding to wind tunnel flow 

ANALYSIS 

Consider the flow over an airfoil spanning the width of a two-dimen- 
sional wind tunnel (figure 1). It is assumed that the boundary-1 ayer thick- 
ness on the sidewalls of the empty tunnel is nearly constant along the 
length of the tunnel and the corresponding value of the displacement thick- 

■ k 

ness is <5 u . With the airfoil present in the tunnel, the boundary-layer 
thickness changes depending on the local pressure field. An estimate of the 
change in the boundary-layer thickness can be obtained from the von Karman 
momentum integral equation 
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( 1 ) 


il + i<Si ( 2 + H .rf) 

dx U dx Plr 

Tne shearing stresses represented by the term on the right hand side of 

equation (1) can be ignored for attached boundary-layers in comparison with 

the pressure gradient effects since the model chord c is much smaller than 

the boundary-1 ayer equivalent length at the model station (References 4, 6 

and 8). Further, since the tunnel walls are normally diverged to compensate 

for the flat plate boundary-layer growth, neglecting the shearing stress 

term would not have an appreciable effect on the results of the following 

analysis. With this assumption and noting 9 = 6 /H, equation (1) can be 

★ 

written in terms of the displacement thickness <5 

2L = £ . (2 * H - M 2 ) ^ (2) 

6* H U . 

The local Mach number M is related to the local velocity U by 

M = U/a . (3) 

Differentiating equation (3) and assuming total temperature to be constant, 
we get 


dU - 1 dM (4) 

U 1+.2M 2 M 

★ 

Combining equations (2) and (4), the variation of 6 with Mach nunber and 
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shape factor can be expressed as 


dS* dH 



2+H-M 2 

1+.2M 2 


dM 

M 


( 5 ) 


The shape factor variation with Mach number for high Reynolds number tran- 
sonic flows can be approximated by (Reference 4) 


H = 1 + (y - 1) M 2 . (6) 

Using this relation, equation (5) can be integrated in closed form to give 
(see Appendix) 

L. = 1 + - 4N|2 

6* 1 + .4M 2 

u « 

* 

<$ u and M^ refer to the undisturbed value of the displacement thickness 

and the free-stream Mach number; usually known from the tunnel test condi- 

★ 

tions. The value of 6 can then be calculated from equation (7) for the 
corresponding local Mach number M. 

Depending on the value of the local Mach number M, the sidewall bound- 
ary-layer displacement thickness will be different from the undisturbed 
★ 

value 5 U . On the upper surface of the airfoil where the flow accelerates, 
the boundary-layer thickness reduces thus increasing the effective width of 
the flow passage and hence changing the local mean velocity. A1 ternatively, 
this flow feature can be viewed as the test airfoil being placed in a 


r l + .2M 2 M ~V 


A + .2M 2 


M 


( 7 ) 
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channel of increased width. Since the disturbances produced by the airfoil 
persist for long distances in the lateral direction at transonic speeds, the 
extent of thinning of the sidewall bound ary- 1 ayer in the region of the air- 
foil can be expected to be nearly the same along the test section height in 
the proximity of the airfoil. Hence the passages above and below the air- 
foil can, to a first approximation be considered as channels with parallel 
side walls. Also, generally the sidewall boundary-layer thickness will be 
much larger than the airfoil boundary-layer except near the trailing edge. 
Therefore, the corner effects of the airfoil/sidewall junction are not 
likely to be significant. The changes in the airfoil pressure distribution 
due to sidewall boundary-layer effects will be largely influenced by the 
change in thickness in the vicinity of the surface. Hence to a first 
approximation, the effective change in the free stream Mach number at a 
representative chordwise location on the airfoil can be obtained by con- 
sidering the conservation of mass flow through a two-dimensional strip 
across the tunnel width. Equating the mass flows in the strip, for condi- 
tions with and without boundary-layer thinning effects, we get 

(b - 26*) H. (b - 2 6*) M c 

= ( 8 ) 

(1 + -ZMg,) 3 (1 + -2M|) 

★ 

where M is the corrected Mach number. With S estimated from equation 
c 

(7), M c can be calculated from equation (8). For small changes in the 
area, the correction AM to the free-stream Mach number can be obtained by 
simplifying equation (8), to give 
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AM = 

H. 


(9) 


★ 




1 + ,2M 2 

oo 

1 - M 2 

1 1 oo 


SMALL DISTURBANCE APPROXIMATIONS 


The variations in the width of the flow passage from the free-stream 
conditions is small for thin sidewall boundary-layers and its effect con- 
sists significantly of introducing perturbations in the stream direction. 
The associated perturbations in the vertical and spanwise directions are 
assumed to be small. This assumption implies that there is a core of uni- 
form flow between the tunnel sidewalls where the variation in the flow quan- 
tities across the tunnel width is quite small and the airfoil aspect ratio 
effect can be ignored. From equation (2), we have 


*. . ® - (2+H-tf)4! 

6* H U . 

Using the approximate variation for the shape parameter given by equation 
(6), the above expression can be expressed in a simplified manner as (Ref- 
erence 4) 


— - - ( 2 + 1/H - M 2 ) f!l! (10) 

5* U • 

* * 
The change d6 in displacement thickness from the undisturbed value 6^ 

can be obtained to a first approximation by substituting for dU, the 

perturbation velocity du due to the airfoil in a two-dimensional flow. 
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( 11 ) 


d5* = - (2 + 1/H - M 2 ) — 6* 

^00 


Let du^ be the perturbation velocity introduced over the two-dimensional 
field in the longitudinal direction because of change in the flow area. 
From one-dimensional isentropic flow relations, u^ is related to the 
change in the flow area dA by 


where 


du b _ 1 dA 


A = b-26 u 


dA = - 2 dS* 


(12) 


Combining equations (11) and (12) and integrating 



U 


CD 


k 

1-M 2 


u_ 

U 

00 


where k is nearly a contant given by 


k = (2 + 1/H - »£) (2 6*/b) . 


(13) 


(14) 


In a wind tunnel flow with sidewall boundary-1 ayers, within the small 
disturbance approximation the longitudinal perturbation velocity u w will 
be the sum of the perturbations due to the airfoil in a two-dimensional flow 
and that due to change in thickness of the sidewall boundary-layer (u^). 
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i .e 


> 


or using equation (13), u can be written as 

u “ u w (1 + — — ) (15) 

1 -*€ • 

The linearized small disturbance equation for the perturbation potential of 
an airfoil in the two-dimensional flow or in a wind tunnel with no sidewall 
boundary-1 ayer effect, is 

(l-M 2 ) * + <J> =0 (16) 

v °° J xx yy v ' 

where is the free-stream Mach number. The airfoil shape and the pres- 

sure coefficient on the airfoil are given by 

y/c = xf(x/c) (17) 

C p = -2u/U oo . (18) 

In a wind tunnel with sidewall boundary-layers, for the same reference 
conditions (U^ and M^), the pressure coefficients on the airfoil will be 
different from the two-dimensional value. Consistent with the small dis- 
turbance approximation, the measured pressure coefficients (Cp w ) are 
related to the perturbation velocity u w by 
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( 19 ) 


Combining equations (15) and (16), and assuming that the effect of the 
sidewall boundary-layers is one dimensional (i.e. = <J> w ), the gov- 

erning differential equation for the perturbation potential in the wind 
tunnel flow is written as 


( 1 - M 2 ) ( 1 + .) <{> + <J> = 0 

^ CO' V ~ ' T W Y Y Y U I \f\J 


1-M 2 


w,xx 


w,yy 


or 

fi - M 2 + k) d> + <t> =o . 

; y w,xx w,yy 

If the free-strean velocity is U and the airfoil thickness 

oo 

ary condition for equation (21) is 


t. 


( 20 ) 

( 21 ) 

the bound- 


(♦ w ,yly=0 ' U » = U„lf(x/C) . (22) 

The boundary condition (22) follows from the assumption that the sidewall 
boundary-layers introduce negligible vertical and spanwise velocities. It 
may be noted that the equation (21) is exactly the same as that derived by 
Barnwell (Reference 4) with the assumption that the spanwise velocity varies 
1 inearly across the width of the test-section. The above derivation demon- 
strates that the Barnwell's assumption of linear variation in the spanwise 
velocity across the width of the tunnel represents to a first order the 
effect of the change in the width of the flow passage due to changes in the 
sidewall boundary-layer thickness. It can be easily shown that the 
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consequence of a gradual change in area is equivalent to a linear variation 
of the spanwise velocity. Further, the present analysis is applicable to 
either narrow or wide tunnels as long as the model chord is long enough for 
the effects to be uniform across the span. Using the coordinate transforma- 
tion 


5 = x and n = y V 1 + k 
equation (21) can be written as 


(23) 


t 1 ♦*.« + * 0 < 24 > 

where M = M / V 1 + k 

e 00 1 

From equation (22), the corresponding boundary condition is 

(* w> „) „ =0 - U. Tf(x/C) / yjl + k . (25) 

From the above boundary condition, it follows that the flow over an airfoil 
in a wind tunnel with sidewall boundary-layers is equivalent to a two-dimen- 
sional flow at Mach nunber M on an airfoil with the free-stream velocity 

e 

or the thickness ratio reduced by a factor of yj 1 + k . 

a) When the free-stream velocity in the equivalent two-dimensional 
flow is the same as U^, the pressure coefficient is given by 

V - 2 (» 5 W U - • <26 > 
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This implies that the pressure coefficient in the wind tunnel flow at test 
Mach number M and free stream velocity U will be the same as that in 

OO CO 

two-dimensional flow at Mach number M„ and free-stream velocity U but 

e °° 

on a thinner profile with thickness ratio i e = t / yj 1 + k. 

b) Alternatively, if the free-stream velocity in the equivalent two- 
dimensional flow is U g = / yj 1 + k , the pressure coefficient will be 

C p, e - - 2 n=0 ’ \ (27) 

This implies that the flow in a wind tunnel at Mach number M^ will be 
equivalent to that in a two-dimensional flow at Mach number M g on the same 
profile with the pressure coefficients increased by a factor yjl + k over 
the wind tunnel value. 

In summary, it follows that for comparison of the wind tunnel data with 
theoretical prediction methods, the appropriate corrected values of the 
various parameters are 


Mach number 
Pressure Coeff. 
Normal force 


M = M / V 1 + k 

Q CO V 

C = Vl + k . c ... . 

p,e 1 p, wind tunnel 

C = V 1 + k . C ... . 

n,e ’ n, wind tunnel 


(28a) 

(28b) 

(28c) 


For transonic flow, if it is assumed that the one-dimensional effects can be 
superimposed, it is appropriate that the local values of the Mach number are 
used instead of the free-stream Mach number in equations (15) and (16). 
Then, the corresponding small disturbance equation for the wind tunnel flow 
i s 
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(1 - M 2 + k) <p 


(29) 


where 


+ <}> =0 
w,xx w,yy 


#(x,y) = fg 8 ± 

U 3x 


Following the same procedure as in the subsonic case, equation (29) can be 
expressed as 


f 1 - M 2 ) <i> + <f> =0 

x,e J 9 w,?S 9 w,nn 


(30) 


where 


M. = M /Vl + k 

Xr ^ 6 


and 



t* + m m 2 !* 

8 « 


(31) 


where M ff is the effective local Mach number. Equation (30) is similar 

A* y \Z 

to equation (24) and, therefore, it follows that in addition to the result 
given by equation 28a, the local Mach numbers are scaled by the factor 
yjl + k . However, at transonic speeds because of the nonlinearity of 
equation (30), the pressure coefficient values cannot be scaled directly in 
proportion to change in thickness. This can be done using von Karman's 
similarity rule as shown in the next section. Since at transonic speeds, 
thickness distribution is more important than the thickness, equation (30) 
and (31) can be used in an approximate manner to provide a scaling for the 
local Mach number. Equation (30) can also be expressed approximately as 


u - wy » w> „ * 4. w>nn - (y ♦ i) (Myuj * w> „ 02) 
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RESULTS AND COMPARISON WITH 
EXPERIMENTAL DATA 

From equation (7), it follows that to a first order the relative thin- 
ning of the sidewall boundary-layer is only a function of the local Mach 
nimber M and the free-stream Mach number M^. Figure 2 shows the expected 
reduction in the displacement thickness under different conditions. It may 
be noted that a significant reduction occurs even for local Mach nimbers 
slightly higher than the free-stream value and with further increase in the 
local Mach number, the rate of decrease of the boundary-layer thickness is 
much slower. Of particular interest is the shaded region shown in Figure 2 
which corresponds to conditions prevailing in testing of supercritical air- 
foils. The relative reduction in the displacement thickness is about fifty 
percent on the upper surface of the airfoil where the local Mach numbers can 
be about 1.2 for free-stream Mach numbers between 0.7 and 0.8. Due to this 
increased width of the flow passage, the flow development over the upper 
surface of the airfoil which is supercritical and is important in determin- 
ing the drag divergence Mach number, will be equivalent to that in a two- 
dimensional flow at a lower free-stream Mach nimber. The extent of the 
correction for the measured free-stream Mach nimber upstream in a wind tun- 
nel is shown in Figure 3. This correction is obtained by solving equation 
(8). For a given free-stream Mach number M , the correction increases 
with an increase in local Mach nimber. This increase is significant for 
thick sidewall boundary-layers and at higher free-stream Mach numbers. With 
relatively thin boundary-1 ayers, about 2 % of the test section width, the 
correction AM is constant over a wide range of local Mach nimbers on the 
airfoil surface. This feature can be explained by the fact that the flow in 
the wind tunnel is confined and the area change is not significant with 
increase in local Mach number on the airfoil. Since the shock position on 
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the airfoil is sensitive to the free-stream Mach number, the perturbations 
introduced because of the change in the thickness of the sidewall boundary- 
layer can be transferred to free-stream conditions. This procedure is in a 
manner similar to the match point technique used by Kemp (Reference 9) to 
correct for the top and bottom wall interference effect using the measured 
boundary data. Correction for the measured pressure coefficients on the 
airfoil surface can be obtained from small disturbance theory by multiplying 
the measured values by the ratio of the tunnel free-stream velocity to the 
corrected freestream velocity. 

It may be noted that the analysis is valid for airfoil model chords of 
the order or greater than the width of the tunnel (c/b > 1) so that tne 
assumption of spanwise uniformity of the flow over the airfoil is satisfied. 
This condition is easily met in practice with the current trend of testing 
long chord airfoil models in adaptive wind tunnels having negligible top and 
bottom wall interference. With short chord or high aspect ratio models 
(c/b « 1), it is likely that the midspan measurement, will be relatively 
free of the end wall effects. For intermediate aspect ratios, a suitable 
correction may be required. 

The results of the small disturbance theory provides much simpler 
relations for the global effects of the sidewall boundary-layer influence. 
In summary, it was shown that the wind tunnel flow at subsonic and transonic 
speeds over an airfoil spanning the width can be considered equivalent to 
another two-dimensional flow at a different effective Mach number (M g = 
+ k) on a thinner airfoil with thickness = (t /*\]l + k) . The 
equivalent two-dimensional flow can be related to a number of other two- 
dimensional flows using familiar similarity rules for high speed flow and 
results of Barnwell and Sewall can be obtained as particular cases. This is 
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schematically represented in Figure 4. 

a. Subsonic flow: Using the Prandtl-Gl auret similarity rule the 

equivalent two-dimensional flow at Mach number M g on airfoil with 

thickness t can be converted to another two-dimensional flow 
e 

corresponding to Mach number M^ and airfoil thickness x. Then 
the pressure coefficient in the new flow is related to the measured 
values by 


c - 

L p, corrected i 

- yj 1 + k 



(33) 


The result obtained in equation (33) corresponds to the corrections 
suggested by Barnwell in Reference 4. 


Alternatively, if we consider the pressure coefficient 

in the new 

flow to 

be the same as 

in the equivalent flow but at 

a different 

corrected Mach number 

i s rel ated to M by 
e 

M on airfoil of thickness x 
c 

, then M 

c 

> 

I 1 -^c- T e = yj 

^ 1 - M 2 . x 
e 


or ^ 

1 — 1 
1 

II 

> 

(l + k • y/l - M| 

(34) 

or 

II 

M 2 - k 

oo 

(35) 
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For M 2 <k, the corrected Mach lumber M is not defined and this 

oo C 

case corresponds to the result of Reference 3 and repeated later in 
Reference 11. This anamoly is due to the fact that the effects of 
compressibility and that of sidewall boundary layers are of oppo- 
site nature. That is, for the pressure coefficient to remain the 
same, the airfoil thickness t in the new flow being larger than 
x g in the equivalent flow, the Mach number has to be reduced. Or 
from equation (34) it follows that for k << 1 and for low sub- 
sonic flows, the appropriate approximation for M £ is M g rather 
than that given by equation (35) or Reference 3. 
b. Transonic flow: The von Karman transonic similarity parameter K 

in the equivalent flow can be written as 

1 - M 2 

K = i (36) 

O e (Y+l) M|] 2/3 

At transonic speeds, it is only possible to compare flows on bodies 
with different thickness ratios t 1 and T2 at different Mach 
numbers M^ and M 2 . In particular, choosing the thickness in the 
new flow to be t, the corresponding Mach number M c can be cal- 
culated from 


1 - M 2 
e 


1 - M 2 

c 

[t (y+1) M 2 ] 2/3 


or substituting for M g and 
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( 38 ) 


1 - M 2 + k 1 - M 2 
00 _ c 

m4/3 ^4/3 

00 C 


The solution of equation (38) for given M^ and k, gives the 

Mach number M in the new flow on airfoil of thickness t. This 
c 

result is the same as that obtained by Sewall in Reference 6. The 
pressure coefficients in the two flow fields are related by 


C p>c [ (r+l) * C F 3 

t 2/3 


C p [(Y+l) M|] 1/3 
— — 


(39) 


C 


p,c 



(40) 


Equation (38) for the corrected Mach number M can be written in a 
convenient form, as 


M 

c 


(l_ M 2)3/4 


(1-M 2 ) 374 (1+k) 1/4 


(41) 


This shows that for small k, the corrected Mach number M is equal to 
M g ( = ^oo 1 + k) and the corresponding expression for the pressure co- 
efficient is 


C * C 
P*c P 


(l+k) 1/3 


(42) 
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In Figure 5, a comparison has been made of the results obtained using 
different forms of the small disturbance theory for calculating the correc- 
tion to the free stream Mach number due to sidewall boundary-1 ayer effects. 
The results of the present approach show a continuing increase in the cor- 
rection from incompressible to transonic speeds. At higher Mach numbers, 
the difference between the present results and that of Sewall (Equation 38) 
are not significant. Equation (41) suggests that at lower Mach numbers M c 

« (l+k) _3/4 M^ and the trend will be the same for both the present and 

* 

Sewall 's results. For small values of 2 6 u /b, all the different approaches 
give nearly the same results at higher Mach numbers. The various approaches 
represent equivalent flows, either of them can be considered valid. In the 
form proposed in the present report (Equations 24 and 30), the equivalent 
flow has been considered without directly invoking the similarity rules. 
This method of representing the effect of the sidewall boundary layer as 
causing changes in the airfoil thickness as well as in the test Mach number 
is valid from low speeds to transonic flows. At low speeds, the thickness 
effect will be dominant, whereas at transonic speeds the Mach nimber effect 
prevails and affects the drag divergence Mach number. It may be noted from 
Figure 5, that the numerical values of the Mach number correction obtained 
from small disturbance approximations agree closely with the more exact 
calculations (Fig. 3) for the changes in the width of the flow passage. 
Perhaps, this explains for the good agreement reported in Reference 7 be- 
tween the measured pressure distribution and the calculations accounting for 
the sidewall boundary layer effects by the method of Reference 6. 

In Figure 6 and 7, correlation of the normal force coefficient and the 
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local Mach number data using the transformation derived in Equations 28a, b 
and c is presented. The experimental data measured with different sidewall 
boundary layers thicknesses in the Langley 6" x 19" transonic tunnel for 
these figures were obtained from Sewall (Ref. 10). It is seen from Figure 
6b that the corrected normal force coefficient C n \l+k and Mach number 
M A/ 1 + k correlate well for 2 6/b = .028 and .07. For the thickest 

oo ^ (J 

* 

boundary layer 2 6 u /b = .10, the correlation is somewhat poor possibly 

because of the higher order effects not accounted for by the small 

disturbance approximation. The validity of the present approach at 

transonic speeds is further demonstrated in Figure 7 by scaling of the local 

Mach numbers given by equation (31). Remarkably good agreement is obtained 

★ 

for the local Mach number distribution for 2 6^/b = .028 and .07, on the 
upper surface of the airfoil where the assumptions made in the present 
analysis are nearly satisfied. Good correlation of the experimental data 
suggests that equation (28) can be used to correct the wind tunnel data for 
the attached sidewall boundary layer effects. 

CONCLUDING REMARKS 

A simple method is presented to estimate the relative thinning of the 
sidewall boundary layer due to airfoil flowfield in two-dimensional wind 
tunnel testing of airfoils. It is shown that the increase in width of the 
flow passage due to thinning of the boundary-1 ayer can be considered as a 
correction for the test Mach number. Further, using the small disturbance 
approximation, it is demonstrated that the effect of sidewall boundary lay- 
ers as causing changes in both the test Mach number and the airfoil thick- 
ness that the latter effect is dominant at low speeds. Based on this 


21 



analysis, the wind tunnel flow over an airfoil with sidewall boundary layers 
is represented by an equivalent two-dimensional flow with the test Mach 
number M^ reduced by a factor 1/^1 + k . The measured pressure coeffi- 
cients are increased by^l + k at subsonic speeds, and to a first approxi- 
mation by (l+k) 1/3 at transonic speeds, where k is a constant calculated 
using the undisturbed sidewall boundary layer thickness and shape parameter. 
Comparison of the results of the present analysis with the available experi- 
mental data on airfoils showed good agreement. 
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APPENDIX 


d6* = dH _ 2 + H - M 2 dM 

6* H 1 + 0.2M 2 M 


and H = 1 + 0.4M 2 


Integrating A.l. 
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log 6 


log H - 3 
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1 - .2M 2 
1 + .2M 2 


dM 
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log 


C 


where C is the constant of integration 


or 


log 6 = log H + log C - 3 


rfi_ _ ,4M 

J l M 1 + .2M 2 


dM 


* ' M 

log 6 = log H + log C - 3 log 


1 + .2M 2 


or 


6* = C. H. 
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The constant of integration is evaluated by using the condition 

★ 

M = M ro the value of the 6 corresponds to the undisturbed value 
value of C is obtained as 
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(A. 3) 
(A. 4) 

(A. 5) 
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1 + .4M 2 
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(A.6) 


or the ratio of local value of displacement thickness to the undisturbed 
value can be expressed as 
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1 + 0.4M 2 
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1 + .2M 2 


* 1 + 0.4M 2 
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W V 1 + •*£/ 


(A. 7) 
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Figure 1. Typical sidewall boundary-layer growth above a lifting 
airfoil In a two-dimensional wind tunnel. 
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Figure 3. Correction for free-stream Mach number due to sidewall 
boundary-layer thinning in the airfoil region. 
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Figure 4 


Schematic representation of small disturbance theory results for 
sidewall boundary-layer effects at subsonic and transonic speeds. 
















Figure 6a. Measured normal force data on a supercritical airfoil in the 

LANGLEY 6" x 19" TUNNEL with different sidewall boundary-layer 
thicknesses (from Sewall , Ref. 10). 
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Figure 7b. Correlation of local Mach number distribution on NLR-1 airfoil 

(Ref. 10) using present small disturbance approximation (Eq. 31) 
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